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ABSTRACT 

Temperature, density and abundance profiles of the hot intracluster medium (ICM) are im- 
portant diagnostics of the complex interactions of gravitational and feedback processes in the 
cores of galaxy clusters. Deprojection of X-ray data by methods such as PROJCT, which are 
model dependent, can produce large and unphysical oscillating temperatur e profiles. Here we 
valida te a deprojection routine. Direct Spectral Deprojection (DSDEPROJ; I Sanders & FabianI 
I2007h . showing that it solves some of the issues inherent to model-dependent deprojection 
routines. DSDEPROJ is a model-independent approach, assuming only spherical symmetry, 
which subtracts projected spectra from each successive annulus to produce a set of depro- 
jected spectra. 
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1 INTRODUCTION 

X-ray observations of many low-redshift galaxy clusters find 
strongly peaked emission in the centre from relatively high density 
gas with a radiative cooling time below ~ 1 Gyr. Without a compen- 
sating source of heat, this gas should cool rapidly causing a drop in 
the central pressure and a subsequ ent slow inwa rd flow of the over- 
lying gas known as a cooling flow jFabiarJl994t) . This cooling flow 
gas should continue to cool to low temperatures. However, recent 
high resolution X-ray spectroscopy was unable to find the emis- 
sion s i gnatures of gas cool i ng below about ~ 1 keV Peterson et al.l 
I2OO3I : feaastra etal1l2004 IPeterson & FabiaiJbOOa) . In addition. 



the cool gas masses implied b y a cooling flow exc eeded the sum 
of the observed gas masses (eg. lEdge & Frave? 2003) and inferred 



star formation rates by an order of magnitude ( Johnstone et al 
19871: iHicks & Mushotzkvll200 j : iRaffertv et alJl2006llO'Dea et al 



200 a) . A heating mechanism is therefore req uired to prevent signif- 



icant a mounts of cooling from occurring (see lMcNamara & NulsenI 
( I2OO7I) for a review). 

High spatial resolution images of cluster cores from the Chan- 
dra X-ray Observatory have revealed huge cavities, shock fronts 
and cold fronts in the ICM. Cold fronts are sharp boundaries be- 
tween regions of different temperature and density likely caused 
by gas sloshing in cluster cores. The X-ray cavities, or bubbles, 
in the hot gas have been produced by the intera ction of jets 
from the central AGN wi th the surrounding ICM (eg. Fabian et al 



2003VForman et al.l20o4lFabian et alj2006l : lMcNamara & Nulser 
2007). These cavities provide a direct and relatively reliable means 
of measuring the energy injected by the central radio source. How- 
ever, whilst AGN heating has been shown to be capable of energet- 
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ically balancing cooling loss es ( Birzan et al. 2004 ; Raffertv et al.l 
'2OO6; Du nn & FabiM] I2OO6I ; iMcNamara e t al. '2006*). the exact 
mechanisms by which the energy is transported and dissipated and 
how the balance is achieved are still unclear. 

Structures such as shock fronts and cold fronts provide tools 
to study the microphysics and transport processes at work in the 
ICM. These discontinuities allow the measurement of the bulk ve- 
locities of the gas in the plane of the sky and the stability of these 
features provide constraints on thermal conduction, the formation 
of magnetic draping layers and the strength of hydrodynamic insta- 
bilities. Detailed studies of shocks and cold fronts in clusters have 
shown that transport processes in the ICM ca n be easily suppressed 
( lEttori & FabiarJl200(]| : IVikhlimn et alj|200ll) . 

In order to study shocks, cavities and cold fronts in the cluster 
core, and to relate the properties of cool core gas to the strength of 
the AGN feedback, we require a reliable method for extracting the 
core gas properties. A spectrum extracted from the cluster centre on 
the plane of the sky corresponds to a summed cross-section with a 
range of spectral components from the core to the cluster outskirts. 
The spectral properties at any point on the plane of the sky are the 
emission-weighted superposition of radiation from all points along 
the line of sight through the cluster. These superimposed contribu- 
tions from the outer cluster layers can be subtracted off by making 
an assumption about the line of sight extent, usually assuming the 
cluster is spherical, and deprojecting the emission to produce de- 
projected radial profiles in temperature, density and metallicity. 

The most commonly used deprojecti on routine is P ROJCT in 
the X-ray spectral fitting package XSPEC (Arnaud 1996). PROJCT 
takes spectra from a series of concentric annuli and fits each one 
with a set of models to account for all the layers of projected emis- 
sion plus a model for the innermost emitting region. However, PRO- 
JCT has been found to be unstable to changes in the radial binning 
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and produce oscillating temperature profile s that do not relat e to 
physical changes in the gas temperature (eg. lFabian et al.ll2006h . In 
section [21 we use simulated clusters to show that multiphase gas 
can be responsible for the unstable PROJCT deprojections. 

In section[3] we describe a spectral deprojection routine, DS- 
DEPROJ and compare the results with PROJCT for simulated clus- 
ters. In section m we apply both DSDEPROJ and PROJCT to a set 
of three nearby galaxy clusters which contain shocks, cavities and 
knots of cool gas in the core. Finally, in section|5] we show that DS- 
DEPROJ produces stable radial profiles for simulated clusters that 
have multiple temperature components, sharp breaks in tempera- 
ture and density or that are extended along the line of sight. 

We used XSPEC version 11.3.2 (Arnaud 1996J for al l spec- 
tral fitting. The MEKAL spectral model dMewe et J]|l985L [l98^ : 
lKaastrall99AlLiedahl et al J 19951) and PROJCT deprojection model 
used here are the defaults available in that version of XSPEC. 
Photoelectric a bsorption was modelled with the PHABS model 
tealucinska-Ch urch & McCammon 1992). Abundances were mea- 
sured assuming the abundance ratios of lAnders & Grevessel l ll989i) . 

We assume Hq = VOkms"' Mpc"', £!,„ = 0.3 and Clf, = 0.7. 
All errors in fit parameters are Id unless otherwise noted. 



2 PROJCT 

PROJCT fits spectra extracted from a series of concentric annuli si- 
multaneously, to account for each projected component. Each pro- 
jected spectrum is fitted by one or more components described 
by a set of parameters, such as temperature and density. Assum- 
ing spherical sym metry to calculate suitable geometric factors 
( iKriss etal.lll983h . the projected sum of the components along line 
of sights are fitted to each spectrum simultaneously. PROJCT reads 
the inner and outer radii of each annulus from FITS header items 
to c alculate the volume of e ach corresponding shell. 

Ijohnstone et al.l j2005l) tested the PROJCT model on synthetic, 
single-temperature cluster spectra and found that the fitted tempera- 
ture and density profiles agreed with the true profiles. However, the 
deprojected temperature profiles generated by PROJCT have been 
found to oscillate between values separated by seve ral times the 
uncertainties on the values (eg. for the Perseus cluster. lpabian et al] 
I2OO6 ). This oscillation can disappear if different sized annuli are 
used, indicating t hat this is not related to physical changes in the 
cluster properties. iFabian et al.l ( l2006h suggest that this instability 
may be caused by multiphase gas or deviation from spherical ge- 
ometry. 



2.1 Simulated Clusters 

In order to thoroughly test our deprojection routines, we produced 
a set of simulated clusters with a range of properties including 
two-temperature components, sharp breaks in density or tempera- 
ture and an elongated geometry. The generic simulated cluster was 
based on a 500 ks observation of the Perseus cluster, with redshift 
z = 0.0183 and Galactic absorption «h = 0. 1 x 10^" cm^^ . The ra- 
dial dependence of the deprojected parameters were given by the 
analytical approximations: 



3.9 X 10^2 



4.05 X 10" 



7x fl + (r/100kpc)^) 

T = ^ keV 

2.3 + (r/100kpc)^ 



(2) 



2 = 0.35 + 0.01391 

V kpc 



-0.000243 f-^ 
V kpc 



1.031 X lO"*" I — I Z0 forr<121kpc (3) 
Vkpcy 



:O.3Z0 forr> 121kpc 



(4) 



l + (f-/80kpc)' 



1 + (r/280kpc)2) 



0.87 



cm^^ (1) 



The radial depen dence of the density was determined by 
IChurazov etakl ilOOA where the second term, which describes 
the distribution on larger scales, was taken from Ijones 81 FormanI 
[l999). The ex pression for the gas te mperature distribution was 
also taken from lChurazov et al] 2003h . An approximate fit to the 
observed abundance profile of the Perseus cluster was used for the 
abundance distribution. 

The FAKEIT command and a PHABS(MEKAL+MEKAL) model 
in XSPEC were used to create artificial spectra for the two- 
temperature galaxy cluster with the parameters set according to the 
supplied equations and then folded through suitable response files. 
The input parameters for the second MEKAL component were also 
determined by equations [T}0 but the temperature and density at 
each radius were halved. The projected spectra were generated by 
stepping through the cluster radius on the sky and, for each radius, 
summing all the artificial spectra along the line of sight. The final 
simulated projected spectra were grouped into bins with a mini- 
mum of 50 counts, allowing the use of statistics. 

This same simulated dataset was used throughout sections |2] 
and[3]to allow a consistent interpretation of the deprojection meth- 
ods. The enhanced metallicity in the second annulus from the cen- 
tre, which for these noisy core spectra results in lower temperature 
and density best-fitting values, and the fluctuations in the profiles 
from 50— 100 kpc are particular to this simulated cluster (eg. Fig- 
ure[3}. Deprojections of other simulated clusters do not also display 
these specific features. 



2.2 Single-temperature Model 

We fitted these two-temperature synthetic spectra with an absorbed 
single-temperature MEKAL thermal model and PROJCT in XSPEC 
11.3.2. The redshift and Galactic hydrogen column density were 
fixed to the simulated cluster values z = 0.0183 and nu = 0.1 x 
10^^ cm"^. The temperature, metallicity and model normalization 
(relating to density) parameters were allowed to vary and the re- 
sulting profiles are shown overlaid on the true profiles in Figure [T] 
The single-temperature PROJCT profile appears to bounce between 
the two temperature components suggesting that PROJCT tends to 
account for one temperature component in one shell and the other 
component in a neighbouring shell. If we assume that the single- 
temperature spectral model is a good fit to the data, a very mislead- 
ing result is produced. 

As an aside, we note that the bounce in the outermost annulus 
of the density profile (Figure[T) is caused by a different effect. The 
emission extends beyond the outermost annulus of each cluster we 
analysed so that the emitting volume associated with this annulus 
is too small. This causes an overestimate of the projection on to 
the next annulus in, so that the fit to that annulus underestimates 
its brightness. However, the effect is minimal for the annuli beyond 
the outer two because the steep surface brightness profile limits the 
amount of projection. 
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Figure 1. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a two-temperature simulated cluster The 
single-temperature PROJCT fit to the data is overlaid on the true profiles. 

2.3 Fixing Parameters in PROJCT 

The oscillation in tlie temperature profile can sometimes be allevi- 
ated by fitting the annuli sequentially from the outside and freezing 
the parameters of components in the outer shells befo re fitting spec- 
tra fr om shells inside them (PROJCTFIXED; see eg. ISanders et al.l 
l2004h . This prevents the poorly modelled spectra near the centre 
from affecting the results in outer annuli. However, this method 
underestimates the uncertainties on the model parameters because 
the outer shell uncertainties are not included. Figure [2] shows that 
freezing the outer shell parameters does remove some of the oscil- 
lation in the outer shells but that this approach produces particularly 
poor fits for the innermost annuli. 

2.4 Two-temperature Model 

A two-temperature MEKAL model can be fitted to the data if there 
are sufficient counts to produce well-constrained parameters. We 
reduced the number of free model parameters in each annulus by ty- 
ing the abundance parameters together, leaving the two temperature 
parameters and the normalizations free. The redshift and Galactic 
hydrogen column density were fixed to the simulated cluster values 
z = 0.0183 and;iH =0.1 x 10^2 cm^^. 

Figure[3]shows that a two-temperature model provides a much 
better fit to the simulated profiles with = 1 .00 compared to Xv = 
4.56 for the single-component model. The two-temperature model 
closely follows the true cluster profiles except for the two innermost 
annuli. The spectra from the inner annuli contain a large quantity of 
projected emission which, when subtracted off, leaves particularly 
noisy deprojected spectra for those shells. 

In summary, we have shown that assuming a single- 
component PROJCT model for a cluster with two spectral compo- 
nents, produces an oscillating temperature profile with values that 
are unrelated to physical changes in the cluster gas. Similar depro- 
jection methods which assume a spectral model for the deprojec- 
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Figure 2. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a two-temperature simulated cluster. The 
standard single-temperature PROJCT fit and the single-temperature PROJCT- 
FIXED fit are overlaid on the true cluster profiles. There are no vertical error 
bars for the PROJCTFIXED model (green line) 
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Figure 3. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a two-temperature simulated cluster with 
lines indicating the two spectral components of the true cluster (solid 
and dashed blue), single-temperature PROJCT fit (black triangles), two- 
temperature PROJCT fit (solid and open green diamonds). 
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tion will suffer from the same problem. To provide a solution to 
this issue, we present and validate a new deprojection routine, DS- 
DEPROJ, which assumes only spherical geometry. 



3 DSDEPROJ 

DSDEPROJ dSanders & FabianlllOOTi) is a model-independent de- 
projection method which, assuming only spherical symmetry, uses 
a purely geometrical procedure to su btract off the p rojected emis- 
sion in a series of shells (similar to Nulsen & Bohrin gerlll995l) . The 
resulting deprojected spectra can then be fitted with single or mul- 
tiple temperature models in XSPEC to produce deprojected profiles 
in temperature, density and metallicity. 

DSDEPROJ takes as inputs: spectra extracted from a series of 
concentric annuli in a sector of the cluster and suitable blank-sky 
backgrounds. DSDEPROJ then performs the following steps: 

(i) Subtract the equivalent background spectrum from each of 
the foreground cluster spectra 

(ii) Start at the spectrum from the outermost annulus: divide the 
count rate in each spectral energy bin by the emitting volume, as- 
suming it was emit ted from a sectio n of a spherical shell, (geo- 
metric factors from iKriss et al.|[T983h to give a spectrum per unit 
volume 

(iii) Scale the spectrum per unit volume from the outermost an- 
nulus by the volume projected onto the neighbouring inner annulus, 
and subtract this from the count rate in each spectral bin of that an- 
nulus 

(iv) Calculate a new count rate per unit volume in each spectral 
bin for this inner shell 

(v) By moving inwards to each successive annulus, we can sub- 
tract off the projected emission of each outer shell from the inner 
shells and produce a set of deprojected spectra 

We used a Monte Carlo technique to calculate the uncertain- 
ties in the count rate of each spectral channel in each spectrum. 
Each of the input foreground cluster and background spectra are 
binned to have 200 counts per spectral channel so that Gaussian 
errors can be assumed. The deprojection process is repeated 6000 
times, each time creating new input foreground cluster and back- 
ground spectra by simulating spectra drawn from Gaussian distri- 
butions based on the initial spectra and their uncertainties. The out- 
put spectra are the median output spectra calculated in this process. 
The 1(7 errors on the count rates in each spectral channel are calcu- 
lated from the 15.85 and 84.15 percentile spectra. 

3.1 Single-temperature Model 

Figure |4] shows the deprojection of the two-temperature simulated 
cluster spectra with DSDEPROJ and the previous result from PRO- 
JCT. DSDEPROJ produces a smooth temperature profile that is the 
average of the two separate components (weighted by emission). 
The errors in the parameters increase towards the cluster centre 
as each subsequent annulus has an increasing amount of projected 
emission which must be subtracted off. This is particularly apparent 
for the central two radial bins. 

3.2 Two-temperature Model 

The deprojected spectra produced by DSDEPROJ were also fitted 
with an absorbed two-temperature MEKAL model to check that the 
two components of the simulated cluster are correctly reproduced 




Radius (kpc) 

Figure 4. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a two-temperature simulated cluster. The 
single-temperature PROJCT and DSDEPROJ results ai'e overlaid on the true 
cluster profiles. 

(Figure |5j- We compared the result with that of PROJCT by leav- 
ing the two temperature parameters and the normalization free and 
tying the abundance parameters for the two components together. 
The redshift and Galactic hydrogen column density were fixed to 
the simulated cluster values z = 0.0183 and nn = 0. 1 x 10^^ cm^^ . 
We found that both DSDEPROJ and PROJCT recovered the simulated 
cluster profiles with a two temperature spectral model. 

A detailed validation of the DSDEPROJ routine for a series of 
simulated clusters is given in section|5] 



4 DEPROJECTION OF A SAMPLE OF GALAXY 
CLUSTERS 

We applied DSDEPROJ and PROJCT to Chandra archive observa- 
tions of three nearby (z < 0.1) galaxy clusters: Perseus, Hydra A 
and Abell 262. These clusters each have deep archive observations 
with total exposure times of 890 ks, 200 ks and 1 lOks, respectively. 
Deep observations of nearby bright clusters yield a large number of 
counts for each cluster which ensures a detailed deprojection with 
a large number of fine bins and a significant detection of multiple 
temperature components. 

4.1 Data Reduction 

The data were analysed using CIAO version 4.0 beta 2 with 
CALDB version 3.4.1 provided by the Chandra X-ray Center 
(CXC). The level 1 event files were reprocessed to apply the latest 
gain and charge transfer inefficiency correction and then filtered 
for bad grades. The improved background screening provided by 
VFAINT mode was applied where available. The background light 
curves of the resulting level 2 events files were filtered to remove 
periods affected by flares using the CIAO script LC.CLEAN. Back- 
ground spectra were extracted from the blank-sky background data 
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Figure 5. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a two-temperature simulated cluster The 
two-temperature PROJCT (black solid and open triangles) and DSDEPROJ 
results (red solid and open circles) are overlaid on the true cluster profiles 
(blue solid and dashed lines). 



sets available from the CXC and cleaned using the same method 
applied to the cluster observations. These background spectra were 
normalized so that the count rate of the source and background ob- 
servations matched in the 9 — 12keV band. 

For each cluster, we identified a suitable sector for deprojec- 
tion and extracted a series of spectra in concentric annuli centred 
on the surface brightness peak. In the Perseu s cluster we se l ected a 
sector containing the weak shock analysed in lGraham et alj j2008l) . 
For Abell 262 we defined a sector which enclosed the half of the 
cluster where the X-ray emission was most extended. Spectra from 
Hydra A were extracted from complete annuli. Individual sources 
were identified using the WAVDETECT algorithm in CIAO. Sources 
were visually confirmed on the X-ray image and excluded from the 
analysis. 

Each cluster sector was divided up into a series of annuli to 
give a minimum of 3000 counts in each deprojected spectrum. 
This criterion provides enough deprojected counts to allow a good 
spectral fit and measurement of the temperature, density and abun- 
dance in each shell. All spectra were analysed in the energy range 
0.5 — 7keV using XSPEC version 11. 3. 2 and grouped with a mini- 
mum of 50 counts per spectral bin. Response and ancillary response 
files were generated for each cluster spectrum, weighted according 
to the number of counts between 0.5 and 7keV. The cluster spec- 
tra were then deprojected with DSDEPROJ and PROJCT in XSPEC. 
In addition, we fitted the annuli in series using the PROJCT model 
with the parameters in the outer annuli fixed to their best fit values 
(see Section^. 

For the single-temperature deprojections, an absorbed plasma 
MEKAL model was fitted to each spectrum in XSPEC. The redshift 
was fixed to the values given in Table [T] and the abs orbing column 
density was fixed to the Galactic values given by iKalberla et al.l 
( l2005t) . The temperature, metallicity and model normalization were 
allowed to vary. 



The PROJCT two-component model fitted to these real clus- 
ters was not readily able to find a minimum. The large number of 
parameters meant that the fit was poorly constrained and unable to 
converge on a solution. We exclude this model from our analysis 
of the cluster sample. The two-temperature DSDEPROJ model also 
produced poorly constrained parameters when fitted to the cluster 
sample. The use of this model was restricted to the Perseus clus- 
ter, for which the best quality data is available. We used a set of 20 
annuli with a minimum of 190,000counts per annulus to produce 
deprojected spectra which could constrain a two-temperature fit in 
each shell. 



4.2 The Perseus cluster 

The data used for this analysis were first presented in the deep study 
of the Perseus cluster (Abell 426) in Fabia n et~al] ^200^ . The total 
good exposure time from the combined Chandra observations is 
890 ks. This deep observation has revealed details of the complex 
interaction between the central AGN and the surrounding ICM. De- 
pressions or cavities in the X-ray correspond to bubbles of rela- 
tivistic plasma that have been inflated by jets from the nucl eus and 
displa ced the surrounding gas ( Boehringer et al. 1993 ; Fabian et al.l 
I2OOOI) . The outer bubbles are presumably from past activity, having 
detached from the nucleus and risen buoyantly outwards through 
the cluster core. Shocks and cool gas around the inner bubbles and 
ripples in the surrounding gas, interpreted as sound waves, also pro- 
vide a challenge for any deprojection routine. 

For our deprojection analysis we focused on a sector of the 
cluster core containing a section of the weak shock that surrounds 
the inner bubbles (Figure |6j Graham et al. 2008). Spectra were 
extracted from a series of 20 radial bins in this s ector spaced 
0. 1 arc min in radius. In an analysis of a similar region, lFabian et aU 
( I2OO6I) found the deprojected temperatures produced by PROJCT to 
be unstable depending on which radial bins were used. They con- 
cluded that this may be due to non-spherical geometry or a multi- 
phas e gas. 

ISanders & Fabiad ( |2007|) found significant variations in the 
model absorption parameter, hh, across the image of the Perseus 
cluster. Most of this variation is likely due to the cluster's prox- 
imity to the Galactic plane. Leaving the nn parameter free for the 
spectral fits in this analysis did not produce any significant differ- 
ences in the radial profiles. 

Figure |7] shows the results of the deprojection routines with 
single-temperature spectral models applied to the sector of the 
Perseus cluster. In agreement with .Fabian et al. (2006.), we found 
that PROJCT produced an unstable deprojected temperature profile 
with temperatures depending on the positioning of the radial bins. 
This oscillation, from --^ 20 — 50 kpc in the top panel of Figure [7] 
is caused by the multiphase gas in the Perseus cluster core. DSDE- 
PROJ produced a smooth temperature profile, revealing a decline in 
temperature from 4keV to below 2keV. PROJCTFIXED produced 
similar deprojected profiles although this method does not correctly 
calculate errors on the parameters. 

PROJCT and DSDEPROJ produced similar stable deprojected 
density profiles (middle panel of Figure I?}. The sharp increase in 
density in the outermost radial bin for all three methods shows that 
there is still a significant amount of projected emission in that bin. 
The outer two bins should therefore be excluded from any quanti- 
tative analysis. The two annuli around 11 kpc show a drop in both 
the temperature and density, corresponding to a drop in the thermal 
pressure of the gas. There is also a jump in the density at 20 kpc 
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Table 1. Chandra observations included in tliis analysis. 



Cluster 


Obs. ID 


Aim Point 


ACIS Mode 


Total Exposure (ks) 


z 




Abell 262 


7921 


ACIS-S 


VFAINT 


110 


0.0166 


0.0567 


Hydra A 


4969, 4970 


ACIS-S 


VFAINT 


200 


0.0549 


0.0468 


Perseus 


3209, 4289, 4946, 4947, 4948, 


ACIS-S 


FAINT 


890 


0.0183 


0.132 




4949,4950, 4951,4952, 














4953,6139,6145,6146 













' Galactic absorption column density (in units of 10~^ cm ^ ) adopted in this paper QCalberlae£aL^200^. 




Figure 6. Unsharp mask image of the Perseus cluster made from the 0.3 — 
7keV band by subtracting an image smoothed with a Gaussian of width 
lOarcsec from one smoothed by 2.5arcsec and dividing by the sum of the 
two images. The annuli used for deprojection have been overlaid in black. 
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Figure 7. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for the Perseus cluster. There are no vertical 
error bars for the PROJCTFIXED model (blue line). 



which is associated wit h the weak shock. T hese features are dis- 
cussed in more detail bv lGraham et alj feOOSi) . 

Although the metallicity is not as well-constrained as the den- 
sity or temperature, the fluctuations in the radial profile were found 
to be stable for different radial binning. The structure in the metal- 
licity profile lik ely relates to real b l obs of high or low metallicity 
in the ICM (see ISanders & FabianI ( l2007h for detailed metallicity 
maps). 

A multi-component fit to the DSDEPROJ deprojected spectra 
was only possible for the high quality spectra from the Perseus 
cluster. Each of the deprojected spectra were fitted with a variable 
number of thermal components, using an F-test to ensure that the 
addition of each one was statistically significant. An F-test proba- 
bility of less than 0.1 was required for an extra model component 
to be added. For each component, the Galactic absorption and red- 
shift were fixed to the values in Table[T]and the abundance parame- 
ters were tied together. The temperature, model normalization and 
abundance parameters were allowed to vary. 

Half of the radial bins required an extra thermal component at 
either a lower temperature of 0.5 — 1 keV or a much higher temper- 
ature component of > lOkeV. The northern edge of the sector in- 
cludes a section of the filamentary structures seen in the soft X-ray 
(Fabian et al. 2002)- These filamen ts have been f ound to correlate 
with detections of cool gas in Ha jConseUce et al.l[200ll) and CO 
( Salome et al .120061) . Each of the deprojected annuli likely contains 
a component of this cool gas but in differing proportions. 

We found that the higher temperature component could be 
equally as w ell fitted by a power -law component with a photon in- 
dex r ~ 1.5. 1 Sanders et all ( 1200411 found evidence for a distributed 
hard emission component surrounding the core of the Perseus clus- 
ter by fitting a high-temperature thermal component. If the origin 
is hot thermal gas, this material could be from a shock caused b 
a mer ger with the nearby high velocity system tSanders & Fabia: 
l2007h . 



4.3 Abell 262 

The inner regions of Abell 262 also host complex structures, in- 
cluding bright knots of emission and a clear cavity to the east of 
the cluster centre wh i ch is coincident with a lobe of radio emission 
dBlanton et Zll2004l) . lBlanton et alj j2004l) found that the knots of 
structure in the core are located in the same region as optical [N 11] 
line emission, suggesting that the gas is cooling down to tempera- 
tures of ~ lO'* K in these regions. 

For a deprojection analysis, we selected a sector of the clus- 
ter that covers the extended region of emission to the south of the 
nucleus and avoids the --^ 5 kpc diameter bubble cavity (Figure [Sj. 
The effect of cavities on the deprojection was examined in detail for 
Hydra A, which has a much more complex system of bubbles. For 
Abell 262, we focused on the effect of the cool knots of emission. 
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Figure 8. Unshaip mask image of Abell 262 made from the 0.3 — 7keV 
band by subtracting an image smoothed with a Gaussian of width lOarcsec 
from one smoothed by 2.5 arcsec and dividing by the sum of the two images. 
The annuli used for deprojection have been overlaid in black. 



The single-temperature deprojection of this sector is shown in 
Figure |9] PROJCT generated an oscillating temperature profile and 
a sharp discontinuity in the density and abundance at 6kpc. This is 
the only cluster in our sample for which PROJCT produced spurious 
results in all three parameter profiles. The discontinuity at 6kpc 
coincides with the edge of the dense core and the jump at --^ 13kpc 
may correspond to a bright knot of emission at the northern edge of 
the sector (Figure[8](. 

The sharp drop in the density parameter was caused by the 
unphysically large value of the metallicity; the spectrum was dom- 
inated by line emission with only a minimal amount of continuum. 
We attempted to solve this problem by placing constraints on the 
metallicity parameter. The PROJCT deprojection was repeated with 
the metallicity parameter in the third annulus, centred on ~ 6kpc, 
fixed to IZq. However, this merely pushed the discontinuity into 
a nearby annulus. Imposing an upper l imit of twice solar on the 
metallicity parameter in each annulus ( I Johnstone et alJlzOOSl) re- 
duced the drop in density in the third annulus, although this was still 
inconsistent with DSDEPROJ and the metallicity parameter reached 
the allowed maximum of twice solar in seven annuli. The con- 
straints on metallicity still produced oscillating temperature pro- 
files. 

DSDEPROJ produced a steady decline in temperature and a 
smooth density profile with a steeper gradient around 6kpc indi- 
cating the edge of the dense core. We tested for the significance of 
a second temperature component in the annulus centred on 6kpc 
by fitting the deprojected spectrum with a two-temperature model. 
Although the density and abundance parameters were poorly con- 
strained in this model, the F-test probability indicated that the ad- 
dition of a second component at a temperature of 0.79^Q QgkeV 
significantly improved the fit. 
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Figure 9. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for Abell 262. There are no vertical error bars 
for the PROJCTFIXED model (blue line). 



together in XSPEC. The central 1.5 arcsec of the cluster, containing 
the AGN, was excluded from the regions analysed (Figure [TOt. 

The large-scale X-ray morphology of this cluster is gener- 
ally smooth, however the core reveals the presence of cavities 
in the am bient cluster gas caused by outbursts f r om the cen- 
tral AGN JivicNamara et~aDl200(]l: lOavid et al.ll200ll: iNulsen etaP 
I2OO2I ; iNulsen et alj|2005l : Iwise et al.ll2007b . These bubbles in the 
ICM are coincident with radio lobes from the central FR type I 
radio source 3C 218 (McNamara et al. 2000). Nulsen et al. (20o3) 
found that cooler gas extends outward from the centre of the clus- 
ter, beyond the cavities, in the direction of the radio source axis. 
There will therefore be significant quantities of cooler gas in each 
deprojected annulus. 

The single-temperature deprojection of Hydra A is shown in 
Figure[TT] The PROJCT temperature profile for this cluster is partic- 
ularly unstable with large oscillations in the temperature parameter 
throughout the core. Although the DSDEPROJ result also appears to 
oscillate, the errors are consistent with a smooth temperature pro- 
file for the inner radial bins. The density profiles show a smooth 
increase towards the centre of the cluster with no features that can 
be readily associated with the large system of cavities in the X- 
ray emission. The effects of this structure were averaged out when 
using complete annuli. The metallicity profile was approximately 
consistent with a constant value of ~ 0.5Z,;,. 

Two sectors of Hydra A, one containing the northern set of 
cavities and the second about a perpendicular axis, were separately 
deprojected but we were unable to significantly detect the inner X- 
ray cavities in the former. 



4.4 Hydra A 

Hydra A (Abell 780) was observed by Chandra for two pointings of 
looks in ACIS-S (ObslD 4969 and 4970). Spectra were extracted 
from complete annuli in both observations, deprojected and fitted 



5 DSDEPROJ VALIDATION 

Multiphase gas in observed clusters is not evenly distributed but 
tends to be concentrated in cluster cores and, in addition, the ICM 
contains sharp features such as shocks and cold fronts. We have 
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Figure 10. Unsharp mask image of Hydra A made from the 0.3 — 7keV 
band by subtracting an image smoothed with a Gaussian of width lOarcsec 
from one smoothed by 2.5 arcsec and dividing by the sum of the two images. 
The annuli used for deprojection have been overlaid in black, excluding the 
inner 9 annuH which ai'e very finely spaced. The innermost black circle 
corresponds to a radius of 26kpc. 
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Figure 11. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for Hydra A. There are no vertical error bars 
for the PROJCTFIXED model (blue line). 



tested DSDEPROJ on a variety of simulated clusters containing ex- 
amples of these features to check that the expected result was re- 
covered. 

DSDEPROJ calculates deprojected profiles that are stable for 
almost any choice of radial binning. Figure [72l shows the DSDE- 
PROJ deprojection of a simulated galaxy cluster under two different 
radial bin distributions. The two DSDEPROJ deprojections produce 
profiles that are consistent and stable. As we show later, instabili- 
ties only occur when a radial bin contains drastically different gas 
properties, such as both sides of shock. For comparison, we show 
the equivalent result under a deprojection with PROJCT in Figure 

In section [3] we use a simulated cluster based on the Perseus 
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Figure 12. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a two-temperature simulated cluster. The 
single-temperature DSDEPROJ results for two different radial bin distribu- 
tions are overlaid on the true cluster profiles. 
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Figure 13. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a two-temperature simulated cluster. The 
single-temperature PROJCT results for two different radial bin distributions 
are overlaid on the true cluster profiles. 



cluster to show that DSDEPROJ produces smooth, stable tempera- 
ture profiles which accurately reproduce the simulated profile. The 
density and temperature profiles for these simulated clusters are 
particularly flat in the radial region considered so we also tested 
DSDEPROJ on simulated clusters with steep power-law profiles in 
temperature and density given by the equations 



© 0000 RAS, MNRAS 000, 000-000 



Direct Spectral Deprojection 9 




Radius (kpc) 

Figure 14. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a single-temperature simulated cluster with 
a steep power-law density profile. 




Radius (kpc) 

Figure 15. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a single-temperature simulated cluster 
which is elongated along the line of sight. 



T = 1.26/-"-'keV 
= 12keV forr > 2Mpc (5) 

/^-0.6\ 

Me = 0.12 cm"^ (6) 

V kpc J 

Z = 0.3Zq (7) 

where r is the cluster radius in units of kpc. Both DSDEPROJ and 
PROJCT correctly reproduce the steep temperature drop and density 
profile (Figure[74l(. 

We tested the assumption of spherical symmetry by applying 
DSDEPROJ to a cluster that was stretched by a third along the line 
of sight (Figure \T^. DSDEPROJ and PROJCT correctly reproduce 
the expected profiles; the central radial bin is poorly constrained 
because the spectrum contains residuals from the incorrect subtrac- 
tion of the outer layers. Assuming incorrect cluster geometry does 
not produce oscillating profiles in PROJCT (or DSDEPROJ), however 
DSDEPROJ could be modified for the deprojection of non-spherical 
systems. 

We also constructed clusters with a sharp drop or rise in tem- 
perature or density at a particular radius. Figure [T6] shows the de- 
projected profiles of a cluster with a temperature break at 42kpc. 
DSDEPROJ and PROJCT both reproduce the breaks in temperature 
without any additional oscillation in the deprojected profile. How- 
ever sharp changes in density can cause additional ringing in the 
density and temperature profiles (dashed lines in Figure [TTt. The 
density parameter in the shell containing the density break cannot 
account for the two distinct values causing an under- or overesti- 
mation of the projection onto the next shell. This problem can be 
alleviated by shifting the radial bins so that the density jump occurs 
close to the edge of a bin (solid lines in Figure [TtIi. 

As a final note, DSDEPROJ assumes that the response of the de- 




Radius (kpc) 

Figure 16. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a single-temperature simulated cluster 
which has a sharp temperature break. 



lector does not vary significantly, which is the case for the Chandra 
AC1S-S3 used in this work. An additional routine which uses the 
ancillary response files to correct for changes in the effective area 
has been incorporated into DSDEPROJ, however this is not required 
here. For the ACIS-S3 detector any variation in effective area has a 
negligible effect on the deprojection (Figure [Tst. 
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Figure 17. Deprojected temperature (top), electron density (centre) and 
metallicity (bottom) profiles for a single-temperature simulated cluster 
which has a sharp density break. The true cluster profile is shown by the 
blue fine; the DSDEPROJ and PROJCT deprojected profiles are shown by red 
circles and black triangles respectively. The open points with dashed lines 
indicate a deprojection where the density jump was positioned in the cen- 
tre of a radial bin; solid points and lines denote the deprojection with the 
density jump shifted to the edge of the radial bin. 
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Figure 18. Deprojected temperature (top), electron density (centre) and 
metalhcity (bottom) profiles for the Perseus cluster. The deprojected pro- 
files were created with (blue line) and without (red line) an ancillary re- 
sponse correction to illustrate that there is no significant difference. 



6 SUMMARY 

We have investigated several issues regarding the reliable depro- 
jection of galaxy clusters and validated the deprojection routine 
DSDEPROJ. Deprojection methods that assume a spectral model 
for the deprojection, such as PROJCT in XSPEC, have previously 
been found to produce temperature profiles which bounce between 
unphysical values separated by several times their uncertainty. We 
have shown that this effect is caused by fitting multiphase gas with 
a single-temperature model. Although a two-temperature PROJCT 
deprojection recovers the true profiles for the simulated clusters, 
the majority of the cluster observations currently available in the 
Chandra archive do not have a sufficient number of counts to con- 
strain the two-components in a detailed deprojection. 

DSDEPROJ, our deprojection routine which assumes only 
spherical geometry, solves some of these issues inherent to model- 
dependent deprojection routines. DSDEPROJ produces a set of 'de- 
projected spectra' which can then be fitted by a suitable spectral 
model in XSPEC. We have shown that this method does not generate 
the oscillating temperature profiles for multi-temperature clusters 
and produces a stable solution for an elongated cluster and clusters 
with breaks in temperature or density. 

PROJCT and DSDEPROJ were applied to a small sample of 
three nearby galaxy clusters, each of which has a component of 
cooler gas or complex substructure in the core. PROJCT produced 
a rapidly oscillating temperature profile for each cluster and, in 
Abell 262, unphysical values of density and abundance. The de- 
projected profiles generated by DSDEPROJ were smoothly varying, 
producing stable solutions at shocks, knots of substructure and with 
different radial binning. DSDEPROJ was only able to constrain a 
two-temperature spectral model for the Perseus cluster deprojec- 
tion. This revealed a low temperature component in several annuli 
that contain cool gas filaments and a power-law co mponent asso- 
ciated with a distributed hard emission component jSanders et al.l 
[2004). 

The DSDEPROJ source code is available at www- 
xray.ast.cam.ac.uk/papers/dsdeproj 
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